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The four-time correlation function of a general dynamical variable obeying Gaussian statistics 
is calculated for the trap model with a Gaussian density of states. It is argued that for energy- 
independent variables this function is reminiscent of the four-time functions that have been dis- 
cussed earlier in the interpretation of the results of four- dimensional NMR experiments on su- 
percooled liquids. Using an approximative relation between the four-time correlation function and 
the cubic response function the nonlinear susceptibility is calculated and the results are compared 
with the corresponding ones resulting from an exact calculation. It is found that the results of the 
approximation change the qualitative behavior of the modulus of the susceptibility. Whereas in the 
exact calculation a peak is found in the modulus in most cases, depending on temperature and the 
additional model parameters no such peak occurs in the approximation. This difference has its 
origin mainly in an incorrect estimate of the static response. The results are discussed in relation 
to recent experimental findings. 

I. Introduction 

In the last decade a number of steps have been undertaken in order to deepen our understanding of 
the heterogeneous dynamics of supercooled liquids and glasses [U EJ El HJ [5] . In many investigations 
certain higher-order time correlation functions play an important role. This holds for computer 
simulation studiesjEJ [7], frequency-selective techniques like nonresonant holeburningjH], [9j E] and 
four-dimensional NMR experiments [Till [TT1 H~2] . Apart from these methods to characterize the 
dynamic heterogeneities it has been shown how a length scale can be extracted from a special four- 
point correlation function x* (t) [IB] and its detailed properties have been studied theoretically [141 
IT5| [TB] . It is to be noted that the special function Xi(f) differs from the corresponding four-time 
correlation function that had been used earlier in the experimental determinations of the length 
scale associated with the heterogeneities of supercooled liquids [T71 ITS] . 

Since it is not straightforward to investigate x*(t) experimentally, it is necessary to have 
observables at hand that are related to it in a clearcut manner. In order to provide a means for 
the experimental determination of the number of correlated particles, N corr , or the length scale, 
Bouchaud and Biroli related the nonlinear (cubic) response Xsi^iT) to a four-point function |19j. 
This relation was utilized in the experimental determination of N nnrr [201 |2T] . It was argued that 
the modulus of the cubic response function, |% 3 (o;, T)|, exhibits a hump-like structure which is 
assumed to be a distinctive feature of glassy correlations [20]. The maximum of |x 3 (o;,T)| is 
expected to decrease with increasing temperature and to be directly proportional to N COIT . In 
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systems without glassy correlations one does not expect a peak but a 'trivial' behavior, i.e. a 
smooth decay to zero starting from a finite low-frequency limit as a function of frequency. 

Recently, nonlinear dielectric experiments have also been performed in a different context 
to investigate the dynamic heterogenenities of supercooled liquids [22], [23] and to measure the 
configurational heat capacity of liquids [24|. 

A nonlinear response theory for Markov processes has been presented in ref . [23J , to be denoted 
as I in the following. The general theory was applied to two simple models for glassy relaxation, 
the model of dipole reorientations in an asymmetric double well potential (ADWP-model) [2B1 127] 
and the well-studied trap model with a Gaussian density of states[2Sl ESI GEH ED E2J E3]- The 
results of the model calculations suggest that a direct relation between the cubic response function 
and some type of glassy correlations cannot be shown to exist in general. For the ADWP-model, 
a peak in the modulus of the cubic response is observed only in a certain narrow temperature 
range and only for non-vanishing asymmetry. For most temperatures the model exhibits 'trivial' 
behavior. The trap model shows both, a peak or trivial behavior, depending on the variable 
chosen and on temperature. Furthermore, for some specific choice of the dynamical variable used 
to probe the dynamics, the peak-maximum increases as a function of temperature and for other 
choices it decreases. Both models are of a mean-field type and therefore of course do not show 
any aspects of spatial correlations as they exist in finite dimensional systems. 

The general theory provides expressions for the nonlinear response functions that cannot be 
related to well-defined time correlation functions in any obvious manner. In the theoretical paper 
by Bouchaud and Biroli[19j, it was argued that for systems obeying a Langevin dynamics, the 
cubic response near a phase transition should be related to the four-time correlation function via 
an expression that reminds of a fluctuation-dissipation theorem (FDT): 

Rs{to,t u t 2 ,t 3 ) ~ (3 \ f (M{t)M{h)M{t 2 )M{h)) (1) 
at\dt 2 at 3 

Here, M(t) is the dynamical variable conjugate to the applied field, t > t 1 > t 2 > t 3 , 
R3(to,ti,t 2 ,t 3 ) is the (impulse) response to field kicks taking place at ti, t 2 and t 3 and (3 = 1/T 
with the Boltzmann constant set to unity. This relation is meant to indicate that the right hand 
side is the dominant contribution to the response and that no other term contributing to it will 
be more divergent. 

In the present paper, I will calculate the four-time correlation function occuring in eq.(JT]), 
investigate its relation to the response and furhermore show that for the special case of a trap model 
this function is directly related to the class of four-time functions observed in four-dimensional 
NMR. In the next Section, this relation will be discussed in more detail and afterwards the cubic 
response originating from eq.Q will be compared to the exact expressions given in ref. [23]. Finally, 
some conclusions will close the paper. 

II. Four-time correlation functions for trap models 

In order to calculate correlation functions of a dynamical variable M(t) obeying a Markovian 
dynamics one uses the relevant probability functions as obtained from the dynamical rules. The 
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stochastic dynamics for the trap model is defined by the master equation (ME) for the conditional 
probability to find the system in the trap characterized by the trap energy e at time t provided it 
was in trap eo at to, G(e, t + toko, *o) — G(e, t|eo, 0) = G(e, t|eo): 

G(e, t\eo) = -K{e)G{e, t\e Q ) + p{e) J de' K (e')G(e' , t |e ) (2) 

In eq.(|2]), the escape rate is given by 

«(e) = K^e* (3) 
with the attempt rate k^. The model with a Gaussian density of states (DOS) is defined by 

P(e) = ^e-^ (4) 



with a = 1. Note that in this model the system reaches equilibrium at all temperatures T 
(measured in units of a) and the equilibrium populations are found to be Gaussian p eq (e) = 
lim^oo G(e, t\e ) = ~^^ e ~^~^' l ^ 2a '^ with e = — (id 1 . In contrast, the often studied model with an 
exponential DOS shows a transition to a low-temperature phase in which no equilibrium is reached 
and the system ages for all times[29]. In the following calculations I will always assume that the 
system is in thermal equilibrium. Thus aging effects are unimportant in the present context. 
The two-time correlation function (2t-CF) of a variable M(t) in general is given by: 

C M (t,t ) = (M(t)M(to)) = Jdej de M(e)M(e )G(e,t- t \e )p cq (e ) (5) 

As in I, a Gausian approximation for the correlations of the dynamical variables M(e) will be 
used, 

(M(e)) = and (M(e)M(e )> = 5(e - e )(M(e) 2 > (6) 

yielding C&f(t, *o) = Jde(M(e) 2 )G(e,t — t |e)p eq (e). For the trap model, the probability that the 
system returns to a trap that has been occupied earlier vanishes. If there are N traps, this means 
that one has G(e,t|e ) = 6(e - e )e~ K{e)t + 0{l/N). One thus finds 

n 2 (t) = j de(M{e) 2 )e- K(e)t p cq {e) (7) 

The most obvious interpretation of this function is obtained for (M(e) 2 ) = 1, in which case II2 (t) 
is directly related to a dynamic structure factor at very large values for the scattering vector [29]. 
This means that every movement completely decorrelates the variable. A similar function also is 
observed in NMR when stimulated echo techniques are used[2] with a high spatial resolution, i.e. 
a large value of the evolution time t p . In that context, n 2 (t) has been termed the angular jump 
function F 2 (t)$M\. 

Next, the four-time correlation function (4t-CF) occuring in eq.(JI]) will be calculated, 

C 4 (t,t u t 2 ,t 3 ) = {M{t)M{t x )M{t 2 )M{h)) (8) 

In order to proceed, I will again use the Gaussian factorization approximation 

(M(e 1 )M(e 2 )M(e 3 )M(e 4 )} = 6(e 1 - e 2 )6(e 3 - e 4 )(M(e 1 ) 2 )(M(e 3 ) 2 ) 

+ 5(6, - e 3 )5(e 2 - e 4 )(M( ei ) 2 )(M(e 2 ) 2 ) (9) 
+ 6(e 1 - e 4 )5(e 2 - e 3 )(M( ei ) 2 )(M(e 2 ) 2 ) 
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For the trap model, this yields 



Ci(t, h, t 2 , t 3 ) = n 4 (t, t u t 2 , t 3 ) + 2^ 2 {t - t 3 ) (io) 

* 2 (t) = / rf e (M(e) 2 ) 2 e- K(e) y q (e) (11) 



where I defined 



and the specific 4t-CF 

U4t ,h,t 2 ,t 3 ) = JdeJde\M(ef)(M(e') 2 )e-^G(e,t 12 \e^ (12) 

In this last expression, I used the abbreviation tki = (t^—ti) to indicate that only the time intervals 
are relevant, 

n^O; ti, t 2 , t 3 ) = 114(^01,^12^23) (13) 

The special dynamic character of the trap model becomes particularly apparent if one uses energy- 
independent variables, (M(e) 2 ) = 1. Setting t i = ^23 = tp and £12 = t req , one recovers the 
4t-function F^f, t req ) as it results in the framework of a free-energy landscape model for glassy 
relaxation, if the assumption of large t p is used, cf. eq.(16) in ref . [35] . 

F A (t F , t req ) = JdeJ de'e-<^G(e, t req \e')e-< e '^p e \e') (14) 

The interpretation of this function is quite straightforward. In the first time interval, one selects 
slow molecules, i.e. only those molecules which have not reoriented during tp (due to the factor 
e -n{t')t F - m e q (14)). i n the 're-equilibration' period, the system evolves freely, i.e. transitions 



among the different traps can take place (factor G(e,t req \e')) and finally one only detects those 
molecules that still are slow (e~ K ^ tF ). The fraction of selected molecules is determined by the 
limit of long t req , denoted by F^tp, 00). Similar arguments can be set up for another prominent 
4t-CF, which has been denoted G^(tF,t req ,t) ~ n 4 (tp, t req , t). This function is not designed to 
detect the 'lifetime' of the dynamic heterogeneities, but the rotational decorrelation of a selected 
subensemble[34j. 

It is well known that the 2t-CF II2 (t) does not obey time-temperature superposition and the 
same holds for the function F 4 , as can be observed in FigJI] In FigjT^i, -F 4 (tp, t req ) is shown for two 
temperatures as a function of the re-equilibration time. Here, r eq = n^e^l 2 ^ a is the integral 
relaxation time of H 2 (t) for n = 0, r eq = /(ieK(e) _1 ]9 cq ( e )- Note that F^f, t req ) decreases more 
exponentially at higher temperature, as shown for tp — 1 in the inset of Figjpi. This is similar 
to the 2t-CF and is a manifestation of the failure of time-temperature superposition in the trap 
model with a Gaussian DOS. FigjljD shows the time scale of the re-equilibration, denoted by r 4 . 
This time scale T4 is an increasing function of the filter efficiency FEj2]. The steepness in this plot 
simply reflects the different width of the underlying relaxation time distribution. All features are 
completely compatible with the fact that a given relaxation rate looses its memory on the time 
scale of the relaxation itself [36J. Of course, for the trap model there is only one time scale for 
relaxation and thus this behavior is to be expected. 

In a next step, consider an energy-dependent variable. As in I, I choose the following form[ 



(M(e) 2 ) = e" n/3£ (15) 
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W T -q FE=1-n 2 (t p ) 

Figure 1: a: (left) F^(tp,t req ) as a function of t req for temperatures T jo = 0.4 and T/a = 1.0. The 
filter times have been chosen as tp = 10~ 2 (full lines), tp = 10 _1 (dashed lines), tp = 1 (red (grey) 
lines), tp = 10 (dotted lines) and tp = 100 (dot-dashed lines). All times are measured in units of 
r eq . The inset in the lower panel shows F±{tp,t req ) — F 4 (£p, oo)/(l — F^(tp, oo)) for tp = 1 and 
T = 0.4cr (dashed line) and T = a (full line), b: (right) T4/r cq as a function of the so-called filter 
efficiency, FE[2], where r 4 denotes the 1/e-decay time of F^(tp,t req ). 



with variable n and where the static value of M 2 has been set to unity. This choice leaves the 
2t-CFs unaffected, because of 

nPe-n{e)t _ „2&£?>0V f , , x««-K(e)t 



dep(e) eq e- npe e- Ki - £)t = 6 / rfep(e) eq e 

with t n = t ■ e~ nl32(T ' 2 . Thus, only the time scale and the amplitude of n 2 and ^ 2 are changed, cf. 
also the discussion of the linear susceptibility in I. 

However, the interpretation of F 4 (tp,t req ) in terms of selected subensembles can be flawed. If 
energy-dependent variables are considered, the factor (M(e') 2 )e~ K ^'- ) ' F determines the dependence 



of the corresponding function H^tp, t req , tp) on the 'selection' period, cf. eq.(12). If for instance, 
n = 1 in eq.(15), then sub-ensembles with large rates (e /3e 3> 1) have very little weight in the 



selection (e _/3e 1). Thus, while for n — 0, the plateau-value F±(tp, oo) is a monotonously 
decaying function of tp, this must not hold in general. In Figj2j this plateau-value is shown for 
T = a and various values for n. It is obvious, that for n < the fraction of selected molecules 
apparently first increases and decreases only for long tp. This fact shows, that the character of 
the 4t-CFs can change dramatically for different variables, even though the 2t-CFs are unaffected 
by the choice of n. 

Note that in the context of the four-dimensional NMR experiments, only n = is relevant 
and there is no ambuigity. In general, however, the 4t-CFs cannot be interpreted in the simple 
picture and also the fact that there is no return probability in the trap model does not allow a 
simple interpretation. This situation is very similar to the findings for the response, where the 
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Figure 2: The plateau value F 4 (t^,oo) as a function of the filter time tp vor various different 
energy-dependencies of the dynamical variable, i.e. different values of n, as indicated. The red 
(grey) line represents n — 0. 

linear response does not change its features with n but the cubic response strongly depends on 
the choice of n. 

III. Nonlinear response functions from 'quasi-FDT' 

As mentioned above, in ref.[19] it was argued that the most important behavior of the cubic 
response function is captured by eq. ([!]). Given the 4t-CF, we are in the position to compare the 
results of this 'quasi-FDT' relation to the exact response for the same model. 

Before discussing the results of model calculations, the most prominent experimental findings 
are summarized as already indicated in the Introduction. One finds that there is a maximum in 
the one-cu (a — 1) or three-cu (a = 3) component of the modulus of Xzi^), 



with Axi denoting the static linear susceptibility (which in dielectric spectroscopy corresponds 
to Ae). This maximum is a monotonously decaying function of temperature and its value has 
been related to the number of correlated particles or a length scale associated with the primary 
relaxation of the system [20| [21] . 

The findings for X3 as obtained in I are summarized as follows. As for a general Markovian 
dynamics that follows a ME, the response for the trap model depends on the details of the coupling 
of the external field to the transition rates in the sense that it is not irrelevant whether the field is 
coupled to the initial or the final state of a transition from one trap to the other. This dependence, 
however, is weak and does not alter the qualitative behavior of X a (u). Much more important is 
the mentioned dependence of the cubic susceptibility on the dynamical variable chosen, i.e. on the 
value of n. For n = 0, i.e. variables that are independent of the trap-energy, one finds a maximum 
in X a (u). In contrast to the experimental findings, however, this maximum X^ 11 ^ = X Q ,(w max ) 
increases with increasing temperature. For n — 1, one finds 'trivial' behavior, i.e. a decay from 
a finite low-frequency value to a vanishing susceptibility at high frequency, for low temperatures. 




a = 1, 3 



(16) 
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At some finite temperature a peak evolves in X a (ui) with a more or less temperature-independent 
height. For n — — 1, one finds a maximum at all temperatures and one observes that A / ^ max ) 
decreases with increasing temperature, reminiscent to what is observed experimentally. 

Now, eq.Q will be used to calculate the cubic response and the results are compared to the 
exact results as obtained in I. According to eqns.d8|) and (10), one has 



R 



3.FDT 



(to, ti, t 2 , t 3 ) — (3" 



dt 1 dt 2 dt : 



-C^to, ti,t 2 , £3) — /3 l 



d 3 



dtidt 2 dt- 



-n 4 (t , ti, t 2 , £3) 



The structure of II4 given in eq.(12) shows that there is a constant long-time level with regard 



to the second time interval t\ 2 - This is because G(e,t\e') tends towards p eq (e) at long times. 
One can thus write G(e, t\e') = p cq (e) + 5G(e,t\e') where 5G(e,t\e') decays to zero for long times. 
This gives rise to the decomposition n 4 (t i, £12, ^23) = n^(t i,t 23 ) + 5n 4 (toi, t± 2 , t 23 ). For the 
response it is sufficient to consider , f , (5n 4 (toi, £12, £23)- The general expression for the cubic 



dtidt 2 dt 3 " xi *V "U±, -iz, 

susceptibility is given in the Appendix and the steady state result for an external field of the form 
H{t) = H cos (out) can be written as: 



X3,FDT(t) 



E 



=1,3 



(17) 



where H is the amplitude of the applied field and the X3fdt( u ) denote the corresponding 
frequency-components. For a direct comparison, in Figj3] the cubic susceptibility is shown for 
n = 0. It is evident, that Xx } fdt(u) shows trivial behavior, whereas Xi(cj) shows a peak. From 
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Figure 3: Cubic susceptibility x^ 1 ^) according the the quasi-FDT relation, eq.(l) (upper panel) 
and according to the exact calculation (lower panel). Full lines represent the imaginary part, 
dashed lines the real part of x*i\u) and the dotted lines are Xi(uj) using eq.(16). 



the imaginary and the real part of the susceptibilities it is apparent that the dispersive behavior 
in both cases is observed on the time scale of the relaxation time r eq . The main difference lies in 
the zero-frequency limit. As an example, in Fig|l| X a ^FDT(oj) is shown for n = 1 and for different 
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Figure 4: X a> FDr{w) as a function of frequency for various temperatures, T/cr = 1, 1.5, 2, 2.5, 3 
from top to bottom. 



temperatures and only 'trivial' behavior is found. This is at variance with X a (u), where one 
observes trivial behavior at low temperatures and a peak at higher temperatures, cf. Figs. 7 of I. 
This trivial behavior of X a>FDT (uj) is found independent of the value of n and of temperature. As 
already mentioned, the low-frequency limit X a FDT (0) is finite in all CclSGS, clS exemplified in Fig|5| 
Note that for n = and for n = —1, the exact calculation gives X a (0) = 0. For n = 1, the finite 
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Fi gure 5: X at FDT{0) as a function of temperatures for different values of n. The red (grey) line is 
for n = 0. For comparison, for n — 1, the exact result also is shosn as dotted line. 

value of X a (0) does give rise to trivial behavior of X a (u) only at low temperature, but due to its 
decrease as a function of temperature, at higher temperatures a peak is observed. 

These results show that when the dependence of the cubic response on the dynamical variable 
is discussed it is important to note that the occurence of a peak in X a (cj,T) strongly depends 
on the low-frequency limit X a (0,T). Thus, while eq.Q gives rise to a response which resembles 
the dynamic behavior of the exact response, for the trap model the results using this 'quasi-FDT' 
relation give the wrong values for static response with the consequence that only trivial behavior 
of X a ^DT{u) is observed. These observation hold for a = 1 and a = 3 and as in the exact 
calculation one has Xi t FDT(0) = 3X 3j ^dt(0) for the low-frequency limit. 
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V. Conclusions 



I have computed the four-time correlation function for a trap model with a Gaussian density 
of states for energy-dependent dynamical variables using a Gaussian factorization approximation. 
Due to the fact that in this model the probability of returning to the original trap vanishes, this 4t- 
CF for energy-independent variables strongly resembles the 4t-function that has been considered in 
many NMR investigations of the life time of the dynamic heterogeneities in supercooled liquids. If, 
however, one allows for an energy-dependence of the dynamical variables, the interpretation of the 
corresonding 4t-CF no longer is straightforward because subensembles that relax at different rates 
have different weight in the average. The specific choice made in eq.(15), (M(e) 2 ) = exp (— n(3e), 
means that for n = 1 slow subensembles have large weight while for n = — 1 the fast subensembles 
dominate. It is important to note that this different weighting does not affect the 2t-CFs. There, 
only the overall time scale and the amplitude depend on the value of n. In the 4t-CF, the different 
weights prevent the characterization of the time intervals as filter interval and re-equilibration 
interval. As for many other properties of the energy-landscape of glasses, the Gaussian trap 
model captures the main features observed in four-dimensional NMR, particularly the fact that 
there is only a single time scale in the relaxation. This means that the 'exchange-time' is on the 
order of the relaxation time of the primary relaxation or the rate memory parameter is close to 
unity. 

When the cubic response is considered, in general there is no obvious relation to the 4t-CF. 
However, it has been argued that a 'quasi-FDT' relation, eq. ([!]), relating the cubic response to a 
triple time-derivative of the 4t-CF captures the important physics [15]. I have used this recipe in 
order to compute the cubic susceptibility and compare it to the exact results obtained in I. One 
finds that the response takes place on the same time scale as in the exact calculation. This means 
that the approximation basically does not change the dispersive behavior of the response. This, 
however, is to be expected because in the trap model there is only a single time scale. Regarding 
the modulus of the cubic response, the exact results differ strongly from the approximative ones 
based on eq.Q. One finds that for all choices of the dynamical variables and independent of 
temperature the quantity X a (u,T) defined in eq.(16) shows 'trivial' behavior, i.e. a monotonous 
decay from a finite low-frequency value to zero at high frequencies. In the exact calculations, for 
all values of n a peak is observed in X a (u,T), albeit not necessarily for all temperatures. Thus, 
for the trap model the dynamic behavior of the cubic response appears to be described reasonably 
by eq.(jl]) but the static susceptibility is wrong. This makes the interpretation very difficult since 
it changes the character of X a (u,T) completely. 

More generally, due to the dissipative nature of the response on the relaxation time scale, 
one expects a strong variation of the susceptibility in the corresponding frequency range. Only 
if the peak resulting in the modulus from this dispersion is not too small compared to the low- 
frequency limit, a hump located at a frequency on the order of the inverse relaxation time should 
develop. This is the reason why in the Debye model no hump is observed due to a rather large 
zero-frequency value of the response [38]. The same argument applies to the ADWP-model for 
dipole reorientations in the case of vanishing asymmetry. For finite asymmetry, in this model, 
a hump is observed in a certain temperature range around the temperature at which X a (0,T) 
vanishes, cf. the discussion in I. The main difference between the ADWP-model and the model of 
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Brownian rotation lies in the fact that in the former model X a (0,T) depends on temperature for 
finite asymmetry while for the latter model it is temperature-independent. Consequently, these 
two models mainly show trivial behavior. The interplay between the static susceptibility and 
the intensity of the dispersion due to the dynamics determines the behavior of X a (u,T) and in 
particular the possible occurence of a hump-like structure. Consequently, every approximation 
made in the computation of the response must yield the correct behavior at least qualitatively. 

In summary, I have shown that the 4t-CFs of the Gaussian trap model show the features of 
dynamic heterogeneities that have been observed in four-dimensional NMR experiments. The use 
of the 4t-CFs in the calculation of the cubic response is problematic for this model because it does 
not give a meaningful approximation to the exact response. It is important to note that this finding 
cannot be generalized to other models in an obvious way and therefore further calculations are 
required. In particular, it is important to perform calculations for models that follow a Langevin 
dynamics, because originally eq.Q was derived for this case. 

Acknowledgment 

I thank Roland Bohmer, Gerald Hinze, and Jeppe Dyre for fruitful discussions. 



Appendix: Cubic 'quasi-FDT' response for the trap model 

The general expression for the cubic response can be written as (cf. the Appendix of I): 

X3,F DT (t)= fdUHiU) f tl dt 2 H(t 2 ) f t2 dt 3 H(t 3 )R^ FDT (t,t 1 ,t 2 ,t 3 ) (A.l) 

J to J to J to 

In order to compute the Fourier-components, I start from the discrete version of the expression 
for n 4 (t i,ti2,^23), eq.Q, 

n 4 (t i, W23) = E( M ')( M fc 2 ) e " Kftt01 ^(^2)e- K;t23 pr 

k,l 

with Gki(t) = G(ek,t\ei), K k = n{e^) and (Mf) = (M(efc) 2 ). The Greens function can be expressed 
in terms of the eigenvalues A m and the eigenvectors Skm of the symmetrized matrix of transition 
rates [39|, 



G kl (t) 



\ ~~ eq ^kmSi m e m (A-2) 



Here, the eigenvalues A m = — T m all are negative (and the T m are relaxation rates) and for 
the m = one has A = and Sko = Jp^- Consequently, this yields the long-time limit, 
Gki{t 00) = pj?. Thus, the relevant function, <5n 4 (t 01 , t 12 , £23)5 can be expressed as a sum 
including the terms with m > 0. Performing the time derivatives and the Fouriertransform yields 
in the steady state 

xiioA") = ^/? 3 E E (Ml)(Ml)s km s lm p m - K k )(r m - Kihi^fpfpf4 a 2i(") (a.3) 

4 fcj m>0 
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with the spectral functions 



3K fc r^K f + o> 2 (8/^ - 2r m (/c fe + /cQ - T 2 J 

r K k v 2 m + 2 Kfc r m K t + 31^ + 2ag(4^ - r m ) 
~ r m (4 + w2 )( K 2 + w2) (r^ + 4^) 

K k T m Ki - o; 2 (2/t fc + 3r m + 6^) 

T m (n k + 3ki) + 2n k Ki - 6a; 2 
+^1 + ^(4 + 9^) 
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